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ABSTRACT 

Polarized foreground emission is a potential contaminant of attempts to measure the fluctuation 
power spectrum of highly redshifted 21 cm Hi emission from the epoch of reionization. Using the 
Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER), we present limits 
on the observed power spectra of all four Stokes parameters in two frequency bands, centered at 126 
MHz (z = 10.3) and 164 MHz (z = 7.66) for a three-month observing campaign of a 32- antenna 


deployment , for which unpolariz ed power spectrum results have been reported at z = 7.7 (Parsons 
et al. 2014) and 7.5 < z < 10.5 (Jacobs et al. 2014). The power spectra in this paper are processed 


in the same way as in those works, and show no definitive detection of polarized power. This non¬ 
detection appears to be largely due to the suppression of polarized power by ionospheric rotation 
measure fluctuations, which strongly affect Stokes Q and U. We are able to show that the net effect 
of pola rize d leakage is a negligi ble contribution at the levels of of the limits reported in [Parsons et al.| 
( 2014|) and Jacobs et al.|(|2014|). 


1. INTRODUCTION 

Polarized emission at meter wavelengths is a poten¬ 
tially problematic foregrounds for experiments seeking 
to measure the 21 cm Epoch of Reionization (EoR) 
power spectrum. Smooth-spectrum sources occupy a 
well-defined “wedge” in the c ylindrical (fc||,fcj_) sp ace of 
the EoR po wer spectrum (e.g., Morales et al. ]MT^|Pober 
et al.| [2M^, but it has been understood tor sonie time 
(e.g., fUeil et al.||2011| [Pen et al.||2009|) that Faraday- 
rotated polarized emission can contaniinate the mea¬ 
surement of unpolarized EoR emission. When mapped 
into the power spectrum, this Faraday-rotated polarized 
emission generates power which mimics the high /cy emis¬ 
sion of the Eo R, scattering power in to the otherwise-clean 
EoR window. Moore et al.j (2013 hereafter M13) simu¬ 
lated the effects ot the low-level forest of weak, polarized 
point sources leaking into the full 3-D EoR power spec¬ 
trum. M13 used the best measurements available at the 
time, notably the existing (unpolarized) poi nt source sur¬ 
veys (Hales et al. 1988 Cohen et al. 2004) a nd the few 
polarization measurements below 200 iVlHz (Pen et al. 
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2009 fBernardi et al. 2009), and found that the level of 
emission could be probleniatic for any experiment with 
^ 1% leakage from Stokes Q —t /. More recent work, 
discussed further below, suggests that the actual level 
of point source p olarization is lower than M13 assumed 
( Asad et al.|2015 ). Polarization leakage may also be pos- 
sible to address with modest calibration requirements in 


future experiments like SKA (Sutinjo et al. 


2015). 


Polarized emission may arise from diffuse (presumably 
Galactic) emission and point sources. Historically, the 
first measurements of the polarization of diffuse emission 
in the sout hern hemisphere at frequenc i es bel ow 1.4 GHz 
are due to Mathewson & Milne (1964 1965), who used 
the Parkes telescope at 408 M Hz and 0.8"" resolution t o 
follow on t he pr evious work of Westerhout et al. (1962). 
Spoelstraj (|1984|) provided a summary ot the situation at 


frequencies between 408 and 1411 MHz, noting the low 
rotation measures of this emission (< 8 rad m“^) with a 
polarization fraction of ^ 35% near 1.4 GHz and notable 
depolarization at lower frequencies. The emission seems 
to be largely due to nearby sources (^450 pc), and fluc¬ 
tuates on scales from 1 - 1C)°. More recent analysis of the 
an gular power spectrum of diffu se, polarized fluctuations 
by La Porta & Burigana (2006) suggest that the ffuctua- 
tions continue to smaller scale s as a p ower law, depending 
on frequency. Bernardi et al. (2013, hereafter B13), who 


conducted a 2400 square degree survey using t he Murchi¬ 


son W idefield Array (MWA) at 189 MHz, and Jelic et al. 
(2014 hereafter Jel4), both reported significaiit amounts 
of weakly polarized emission on angular scales between 
a few degrees for B13 and a few tens of arcminutes for 
Jel4. The rotation measure (RM) |$| < 25 rad m“^ in 
all cases. 

Regarding point sources, the number of bright sources 
below 200 MHz reported in the literature is small. 
Of particular interest in determining the polarization 
contamination is the linearly polarized fraction p = 


































































2 


Moore, et al. 


^ (Q2 + U'^)/P of sources (where /, Q, and U are the 
Stokes parameters), and, in the absence of low fre¬ 
quency measurements, how this scales from higher fre¬ 
quencies. B13 detected one polarized point source at 
189 MHz, PMN J0351-2744, with a polarization fraction 
p = 0.02 and polarized flux of 320 mJy. At somewhat 
higher frequencies (350 MHz), |Gie 6 ubel et al.| (|2013l 
studied sources behind M31 and found that P350 mhz = 
0.14 P 1.4 GHz, with typical values of a few percent at 350 
MHz. A study of t he depolarization of point sources by 
Fames et al.|(2014|) showed a systematic trend for depo- 
larization of steep-spectrum point sources as frequency 
decreased, resulting in very low polarization fractions 
{p 0.01) below 300 MHz. In a measu rement of M51 
using LOFAR between 115 and 175 MHz, Mulcahy et al. 
(2014) detected 6 background sources with Stokes 1 be- 


tween 44 and 1500 mJy. They detect two “sources” with 
relatively large polarized fractions (0.026 and 0.029), but 
these are the lobes of radio galaxies with unpolarized 
cores. If one were to include the flux from the core, thus 
giving the total polarized fraction integrated over the 
galaxy, this would be lower. For the unresolved sources, 
the polarized fractions are all < 2.8 x 10“^. Depolariza¬ 
tion from 1.4 GHz range d between 0.03 and 0 .2 for the 
six sources. Importantly, Mulcahy et al. (2014) had very 
high angular resolution (20-25") relative to H13 and PA¬ 
PER, thus minimizing beam depolarization effects from 
integrating over sources. 

It is worth noting that the reported rotation measures 
for point sources vary over the range of values re ported in 
the all-sky RM map of Oppermann et al.| (2012), whereas 
the diffuse emission typically shows low RiVl. 'I'his is con¬ 
sistent with the interpretation of this emissio n as nearby 
diffuse Galactic synchrotron. As discussed in |Jelic et al.| 
( 2010 ), polarized synchrotron generated from within a 
magnetized, ionized plasma will be depolarized, and will 
further show structure in the rotation measure spectrum 
of the source. The result is a lower apparent Faraday 
depth. By contrast, extragalactic point sources see a 
large Faraday column through the Galaxy, resulting in a 
high Faraday depth. 

The goal of this paper is to directly quantify the ef¬ 
fect of polarized leakage on upper limits on the EoR 
power spectrum previously from the Donald C. Backer 
Pr ecision Array for Pro bing the Epoch of Reionization 


Parsons et al. (|2014[ henceforth referred to as P14) 


and Jacobs et al.| ( |2U14[ henceforth referred to as J14). 
We begin in bection by reviewing key points of the 
effect of polarization on the EoR power spectrum. We 
describe the dataset and processing steps in Section 
and present power spectra from two frequency bands in 
all four Stokes parameters in Section along with impli¬ 
cations of these measurements for the level of polarized 
leakage into the 21cm EoR power spectrum. We conclude 
in Section [5l 


2 . THE PROBLEM OE POLARIZATION 


In this section, we will briefly review the problem of 
polarization for EoR power spectrum measurements. We 
recall that the index of refraction of an ionized, magne¬ 
tized plasma is birefringent. The left- and right-circular 
polarizations of an electromagnetic wave passing through 
such a plasma undergo different phase shifts, known as 
Faraday rotation, such that the phase difference of 


the light becomes 

= i ^ f neis)Bi\{s)ds = A 2 $, ( 1 ) 

where is the electron density of the plasma, is 
the component of the magnetic field along the line of 
sight, e and rrie are the electron charge and mass, 
is the wavelength of the incident light, and the integral 
extends along the line of sight. Equation defines the 
rotation measure 4). Faraday rotation affects the linear 
components of the Stokes parameters such that a polar¬ 
ized source with intrinsic Stokes Q and 17, when viewed 
through a magnetized plasma, will have measured Stokes 
parameters 

{Q + iU)meas = ^ ^ {Q + iU)intr ( 2 ) 


The frequency structure induced by Faraday rotation 
differs from normal smooth-spectrum foregrounds, and 
exhibits high covariance with the high line-of-sight fey 
modes, which are typically free of synchrotron foreground 


emissio n (e.g. 


Morales et al. 2006 2012 Parsons et al. 


2012b). In fact, there is a nearly one-to-one correspon- 


dence between a 4> mode and the k\\ mode it most infects 
given by 

H{z) 


fell = 4 


c(l -I- z) 


4>A^ 


( 3 ) 


where z is the redshift of the observation, H{z) is the 
Hubble parameter at that redshift, A is the wavelength 
of the observation, and $ is the rotation measure in 
questionj^ At 164 MHz, the central frequency of one 
of the bands we present, a typical rotation measure of 
20 rad m“^ will infect values of around 0.25 /iMpc“^, 
well outside the horizon limit for smooth-spectrum fore¬ 
grounds on short baselines (fc « 0.05 /iMpc”^ for z = 7.7 
and a 30 m baseline). 

The power spectrum for 21cm EoR measurements is 
unpolarized, so the frequency structure induced by Fara¬ 
day rotation must leak into I measurements through in¬ 
strumental effects. We note that any instrument effect 
which leaks Stokes Q or U into / is subject to the kind 
of contamination we discuss here. For PAPER, the par¬ 
ticular form of the dominant leakage comes about as 
follows. Since PAPER has little imaging capability in 
its maximum-redundancy configuration, we cannot form 
Stokes parameters in the image plane, but rather, we 
combine visibilities as if they were images, by defining 
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where Vj,: 

ibilities for each frequency and time, and V/, etc., are 
the “Stokes visibilities”. The detailed expression for the 
leakage of Stokes Q, U, and V into I due to wide-held 
beam s has been explored elsewhere (e.g., Shaw et al. 
2015), and we defer a detailed study of this effect tor 


A similar e quation to this h as been presented in two other 
papers, M13 and|Pen et al.| (|2009|). Both of these contain algebra 
mistakes, which are corrected m the formulation we present here. 
We thank Gianni Bernardi for pointing out these mistakes. 
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PAPER to future work. Considerable simplification oc¬ 
curs in the limit that the two feeds produce orthogonal 
response to the two electric field polarizations everywhere 
on the sp here, which is not strict ly achievable in practice 
(see, e.g.,|Carozzi & Woan|2009|, but a good approxima¬ 
tion over most ot the 45"" EWHM PAPER field-of-view. 
In this limit, expanding the first row of the matrix and 
expressing the visibilities in terms of the beams Ap for 
polarization p, the baseline b, the unit vector s, and the 
intrinsic polarized signals I and Q, we find 


Vi = 

-f J (A,, - Ayy)Qe-^^^^-^’'^^ dfl. (5) 


The expression for Vq is similar, with I and Q inter¬ 
changed. Thus, if the primary beam of each element is 
not symmetric under rotations of 90°, then the I visi¬ 
bility will have a contribution from both / and Q. This 
provides the mechanism for the spectrally non-smooth 
Faraday-rotated polarized emission to infect the typically 
unpolarized 21cm EoR power spectrum for PAPER. 


3. DATA PROCESSING 

The data used to create these results is nearly identical 
to that presented in this paper’s sister papers, P14 and 
J14. We will review the processing steps presented there, 
highlighting slight modifications. 

paper’s 32 antennae were arranged into a 4 row x 8 
column grid during the 2011-2012 season (PSA32), when 
this data was taken. The East-West row spacing was 
30 m, and the North-South column spacing was 4 m. 
This choice of antenna configuration maximizes base- 
line redun dancy, achieving heightened sensitivity (jPar- 
sons et al. 2012 a), and allows for the r edundant calibra¬ 
tion of Visibilities (e.g.|Liu et al.||2010D. 

This data set consists ot data taken continuously from 
Julian date 2455903 until 2455985, for a total of 82 days 
of observation. The effective integration time for any 
point on the sky is shown in Figure [T] Data is not con¬ 
sidered when the sun is above -5° in altitude. 


3.1. Initial Processing 

We begin with an excision of radio frequency interfer¬ 
ence (RFI), a three step process. First, we flag known 
frequency channels containing nearly constant RFI, for 
example, the 137 MHz frequency bin that contains the 
continuous signal from a constellation of communications 
satellites. Away from these regions, the fraction of data 
flagged is ^ 1% of the total. Next, we subtract adjacent 
time and frequency channels from each other to cancel 
the bulk of the signal, and flag 6tT outliers in the differ¬ 
enced data. Finally, we remove a foreground model and 
flag 4(7 outliers in the residual spectra. 

Once RFI excision has taken place, we low-pass filter 
and decimate the data in time and frequency, in the man¬ 
ner described in Appendix A of P14. This reduces the 
data volume by roughly a factor of forty — from 1024 to 
203 channels, and from integration times of 10 seconds to 
34 seconds — and preserves all celestial signal, including 
the EoR. 

Next, we derive a fiducial calibration solution for a sin¬ 
gle day’s worth of data. We begin by solving for antenna- 



1 


Fig. 1. — Effective integration time per pointing, as a function 
of position on the sphere. This metric is (iefined to give the total 
integration time when integrated over position on the sphere. The 
total field of view surveyed is 2.39 sr. RA zero is at the center 
of the Mollweide projection, and increases to the right, wrapping 
around the sphere. The maximum occurs in RA at approximately 
6h, and in Dec at approximately -30°, above the latitude of the 
array. 


based gains and electrical delays which enforce redundant 
measurements across redundant baselines in the array. 
This reduces the calibration solutions to a single, over¬ 
all gain and delay per polarization of the array, which 
we solve by fitting visibilities to a model of Pictor A, 
corre cting for the p r imary beam gain towards Pictor A, 


as m Jacobs et al. (2013). For sources at zenith, cali¬ 
bration errors are a tew percent. We apply this fiducial 
calibration solution to all data. 

We develop a model of the smooth spectrum fore¬ 
grounds for each integration and each baseline by con¬ 
structing a spectrum of delay components over the full 
available bandwidth (100-200 MHz) using a 1-D CLEAN 
algorithm. We constraint these CLEAN components to 
lie within the entire horizon-to-horizon range with an ad¬ 
ditional extent of 15 ns beyond either horizon. This pro¬ 
cedure does not affect high-delay signal due to the EoR, 
but both removes foreground signal and deconvolves by 
an uneven RFI flagging kernel. We subtract this model 
from the data at each integration and baseline. This 
procedure is described in P14 and J14. 

Finally, we remove cross-talk, defined as an offset in 
visibilities constant with time. For each baseline, and 
for each day, we subtract the nightly average of the data, 
enforcing mean zero visibilities. 

3.2. Polarization Calibration 

To begin a discussion of polarization calibration, we 
summarize the redundant calibration procedure we take. 
First, we take the set of xx and yy visibilities and treat 
them like independent arrays. Within each of these ar¬ 
rays, we solve for the Nant antenna-based gains and Nant 
antenna-based electrical delays which force all baselines 
of a certain type to be redundant with a fiducial base¬ 
line in each type. This leaves three calibration terms 
per polarization to be solved for: an overall flux scale, 
and two delays which set the three baseline types to the 
same phase reference. We solve for these three b y fitting 
the redundant visibilities to a model of Pictor A (Jacobs 


et al. 


So 


2013). 


ar, we have made the reasonable assumption that 
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Center frequency 
z/c [MHz] 

Total bandwidth 
At/ [MHz] 

Channel width 
[kHz] 

Redshift 

z 

Collecting area 
Aeff [m2] 

System temperature 
Tsys [K] 

Leakage FoM 
A-/A+ (Eq.jlgl 

Band I 

126 

7.9 

493 

10.3 

4.47 

836 

3.3 X 10-2 

Band II 

164 

9.4 

493 

7.66 

5.80 

505 

2.2 X 10-2 


TABLE 1 

Observational Parameters 


/ dominates the xx and yy visibilities. If we also assume 
the gains and delays are truly antenna-based, then we 
can apply the calibration solutions from those visibilities 
to the xy and yx visibilities. This procedure omits one 
more calibration term which sets the x and y solutions to 
the same reference phase: the delay between the x and 
y solutions, r^y. 

To solve for this delay, we minimize the quantity 

= Y1 \^*7 “ exp{-2TriT^yiy}\^ , ( 6 ) 


where the sum runs over antenna pairs i,j, finding the 
electrical delay which minimizes Vy in the least squares 
sense. This potentially nulls some signal in VV, but we 
do not expect any significant signal in V at these fre¬ 
quencies. This method o f polarization calibration is sim¬ 
ilar to that presented in Cotton (2012), but rather than 
maximizing Vj/, we are minimizing Vy- By assuming 
that gains are antenna-based and that the flux in V^x 
and Vyy is dominated by unpolarized emission, we need 
not correct for gain differences between Vxy and Vyx (the 
relevant gains having already come from the xx and yy 
solutions). 


3.3. Averaging Multiple Days 

As a final excision of spurious signals (most likely due 
to RFI), on each day, we flag outlying measurements in 
each bin of LS T. W e use the measurement of Tgys out¬ 
lined in Section |3.5| to estimate the variance of each bin, 
and flag 3tT outliers. 

If the data followed a normal distribution, consistent 
with pure, thermal noise, then this procedure would 
flag around one measurement in each frequency/LST 
bin, causing a slight miscalculation of statistics post¬ 
flagging. To counteract this effect, we calculate the ratio 
of the variance of a normal distribution, truncated at 
±3cr (97.3%). We increase all errors in the power spec¬ 
trum by a factor of 1.03 « 1/97.3% to account for this 
effect. 

We compute the mean of the RFI-removed data for 
each bin of LST and frequency, creating a dataset com¬ 
prised of the average over all observations for each LST 
bin, literally an average day. We continue analysis on 
this data. 


Earth’s rotation, and A is the latitude of the observa¬ 
tion. We filter the data in time using a boxcar filter 
in fringe-rate space, defined as one on 0 < / < hE<^® 
and zero elsewhere. While this does null some celestial 
emission (roughly the area between the south celestial 
pole and the southern horizon), its effect is small, since 
the primary beam heavily attenuates these areas of the 
sky. We null these fringe rates as an additional step of 
cross-talk removal. 

Finally, we combine the linearly polarized visibilities 
into Stokes visibilities, as in Equation]^ 


3.5. System Temperature 

Once initial preprocessing has been completed, we take 
advantage of nightly redundancy as a final check on the 
data. Since PAPER is a transit array, measurements 
taken at the same LST on different nights should be to¬ 
tally redundant. This redundancy allows us to measure 
the system temperature via fluctuations in signal in the 
same LST bin from day to day. 

Eirst, we compute the variance in each frequency and 
LST bin over all nights of data <Jjy{i/, t), and convert this 
variance into a measure of the system temperature Tsys- 
This measurement is totally independent of the follow¬ 
ing power spectral analysis, and can be used to quan¬ 
tify the level of systematic and statistical uncertainty 
in the power spectra. It complements measurements of 
Tsys from the uncertainties in power spectra in P14 and 
J14. The variance computed in each frequency/LST bin 
cr^jy{r/,t) is converted into a system temperature in the 
usual fashion: 


^sys iy 1 i') 


■^eff ^Jy 


(7) 


where Ag// is the effective area of the antenna, is the 
Boltzmann constant, 5v is the channel width, and tint is 
the effective integration time of the LST bin. 

Eigure shows the measured system temperature for 
each frequency and LST bin collected during the PSA32 
observing season. To further summarize our data’s vari¬ 
ance, we average Tsys{v, t) over the time- and frequency- 
axes. The frequency-averaged system temperature at 
center frequency may be computed as 


3.4. Final Processing 

After visibilities are averaged in LST, a final round 
of cross-talk removal is performed. Again, we simply 
subtract the time average across LST from the data. 

In the penultimate processing ste p, we pass the data 
throug h a low-pass filter in time. [Parsons fc Bacl^ 
([2009 ) and Appendix A of P14 describe the celes¬ 
tial limits of the fringe rate / for drift-scan arrays as 
cos A < / < bEUJi^, where bE is the east-west com¬ 
ponent of the baseline, a ;0 is the angular velocity of the 


trr. ^ - Iau ^’^c)Tsys {p t) dv 


( 8 ) 


where W (zz; is the window function in frequency used 
and the integral span the bandwidth Az/. Eor ou r anal¬ 
ysis, we use a Blackman-Harris window function (Harris 


1978), chosen to maximally suppress sidelobes. A sim¬ 
ilar expression may be written for the time-axis, where 
our window function is simply the number of redundant 
samples in each frequency channel. 
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Fig. 2. — System temperature in Kelvin as a function of lo¬ 
cal sidereal time (LST) and frequency calculated by Equation 
Solid black boxes enclose the area used to compute the power 
spectrum. The persistent RFI mentioned in Section [3.1| as well as 
the band edges, appear white in this figure. 
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Fig. 3.— (Top Panel) Band-averaged system temperature (Equa¬ 
tion!^ as a function of local sidereal time for band I and II in black 
and mue, respectively. The shaded gray region indicates the range 
in local sidereal time used to compute the power spectra. (Bot¬ 
tom panel) Time-averaged system temperature, averaged over local 
sidereal times IhOOm to ShOOm. The shaded gray and blue regions 
show the window functions for Ban ds I and II, respectively. The 
persistent RFI mentioned in Section [3.1| as well as the band edges, 
appear as zeros in this estimate of ot the system temperature. 


for visibility V(iy) and delay mode r, becomes an estima¬ 
tor for T(k), the Fourier-transformed brightness temper¬ 
ature field. The power spectrum may be computed from 
the delay-transformed visibilities via the equation 


p(k) = 


a 2 y x^Y 

2kB ) ilAv 



( 10 ) 


Here, -P(k) is the three-dimensional power spectrum of 
21cm emission, \^/2kB is the conversion from Jy to K, 
H is the solid angle subtended by the primary beam, Av 
is the bandwidth of the observation, X^Y is the factor 
converting cosmological volume in li“^Mpc^ to observed 


volume QAv (taken from Equations 3 and 4 in Furlan- 


etto et al. 
visibility. 


(2006)), and V{t) is the delay transformed 
ihe k modes are determined by the baseline 


vector and the t mode. 

A subsequent covariance removal, described in detail 
in Appendix C of P14, projects the delay transformed 
visibilities into a basis in which the covariance between 
two redundant baselines is diagonal, and then computes 
the power spectrum from the projected delay spectra. 
This procedure produces an estimate of the power spec¬ 
trum for each LST bin and baseline type. To measure 
the uncertainties in the time-dependent power spectra, 
we bootstrap over both redundant baselines and LST 
samples. 


3.7. Accounting for Ionospheric Effects 

In P14 and J14, we ignored the effects of the iono¬ 
sphere, since for Stokes I these are largely changes in 
sou rce position induc ed by refraction (for a recent study, 
see Loi et al. 2015) and these are negligible on the 


large angular scales considered (Vedantham & Koop- 


mans 


2015a|b ). However, both spatial and temporal fluc¬ 


tuations in the Faraday depth of the Earth’s ionosphere 
will have a strong effect on polarized signal. As the total 
electron content (TEC) varies, it modulates the incoming 
polarized signal by some Faraday depth which is a func¬ 
tion of both the TEC for that time, and the strength of 
the Earth’s magnetic field. Though we assume visibilities 
are redundant in LST for the purposes of averaging to 
form the power spectrum, they do in fact have variations 
due to the variable TEC of the ionosphere. Thus, aver¬ 
aging in LST could result in attenuation of signal. We 
are not able to directly image each day and calculate the 
effects of ionosphere variations based on the properties of 
celestial sources, but we are able to estimate the size of 
the effect on the visibilities used in the power spectrum 
analysis. 


Figure shows the system temperature averaged over 
frequency and LST ranges used to compute the power 
spectra. Frequency- and time-averaged system tempera¬ 
tures for the band centered at 124 MHz (hereafter Band 
I) and 164 MHz (Band H) are reported in Table 


3.6. Power Spectra 

We compute t he power spectrum usi ng the delay spec¬ 
trum approach (Parsons et al. 2012b). For short base¬ 
lines, the delay transform of a visibility, defined as 


V{t) = / V{v)e 




(9) 


3.7.1. Characterizing ionospheric behavior over the 
observing season 

Ionosphere data is provided by a number of sources. 
We have used that from the Center for Orbit Determina¬ 
tion in Europe (CODE), whose global ionosphere maps 
(GIMs) are available in lONosphere map EXchange for¬ 
mat (lONEX) via anonymous ftp. lONEX files from 
CODE are derived from ~200 GPS sites of the Interna¬ 
tional Global Navigation Satellite System Service (IGS) 
and other institutions. The time resolution of GODE 
data is 2 hours, and the vertical TEG values are gridded 
into pixels 5° across in longitude, and 2.5° in latitude. 
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SAST 00:00 2012-02-13 



Fig. 4.— Ionospheric RMs over the PSA32 site at local midnight 
on 2012 February 13, during the observing season described here. 
The plot is shown in local horizon coordinates, with zenith at the 
center, and the horizon around the edges. Note the large, smooth 
variation over the FoV. All of these RMs affect the measured visibil¬ 
ities, making accurate polarization calibration extremely difficult. 



RM (rad m ^) 

Fig. 5.— From top to bottom, a histogram of the zenith iono¬ 
spheric RMs over the season, for the transit of LSTs Ih, 4h and 
8h. The large variance at any given LST correspond s to a large 
attenuation of polarized signal, as discussed in Section|3.7.2| 


We have used the core of th e code provided from the 
ionFR packa^f^^ described in |Sotomayor-Beltran et al. 

[2M5t to access the lONEX hies. As written, 
ionFR provides the RM towards a given RA as a function 
of UT at a given latitude (with the CODE two-hour time 
resolution interpolated into hourly values). 

For each day of the PSA32 campaign, we obtained the 
relevant lONEX file. For a fixed LST, we found the UT 
corresponding to that LST at transit for the PSA32 site 
for each day of observation. Using the data in an lONEX 
file and ionFR we were able to calculate the vertical TEC 
and geomagnetic field (and hence RM) values over the en¬ 
tire hemisphere observed by PAPER. The result is a map 
<f>(U) giving the ionospheric RM induced for sources 
in any direction; an example is shown in Figure 

To give an indication of the time variability of the RM, 
as coad ded into the power spectrum for that LST (recall 
Section 3.3), we calculated for each night of observation 
the zenith RM when 3 LSTs (Ih, 4h, and 8h) were at 
transit. These are shown in Figure!^ Over the season, 
there is a large spread of ionosphericRMs for each LST. 
As the relevant phase shift of the Faraday rotated spec¬ 
trum is this clearly varies by more than a radian 

over the range of days coadded. This will lead to a large 
attenuation of polarized power during LST-binning; we 
calculate this attenuation in the next section. Also note 
that there is a decrease in the average value of the RM 
as LST increases. This is expected, given the strong 
corr elation between the day / night cylce an d TEC val¬ 
ues (Radicella & Zhang 1995 Tariku 2015), and given 
that for this observing season, LS'i'=ih corresponds to 
observing times shortly after sunset, while LST=8h is 
always well into the night. 


3.7.2. Calculating ionospheric attenuation of polarized 

signal 

The effect of the ionosphere requires a modification of 
equation!^ to account for t he effect of Equa. t ion[^ which 
following the formalism of Hamaker et al. ([T996|), is of 


http://sourceforge.net/projects/ionfarrot/ 


the form 

U/ = J dne-^^^'^-^’'/^[{A,, + Ayy)I+ 

{Aj;;^ - Ayy){Q cos 4)+U sin (f)] (11) 

where (f) = 2$(U)A^ and $(U) represents the spatial dis¬ 
tribution of ionospheric rotation measures at the time of 
observation. As we have seen, d)(n) is a slowly varying 
function over the PAPER beam. 

If we assume 4>(U) is spatially constant. Equation 11 
can be re-written 


u;= J dne-^”^-^’'/^{A,, + Ayy)i 
+ cos(j) j dCle~‘^'^'-^'‘"'^‘'{Axx - Ayy)Q 
-Fsin(/> J dne-^^^^-^''/%Axx - Ayy)U 

= Vi + cos cj)VQ + sin (j)Vu (12) 

and we can write the LST-averaged visibility as the Fara¬ 
day rotation-weighted sum of otherwise redundant visi¬ 
bilities: 

V = +cos (j),VQ+sin (fiVu, (13) 


where (ft is the zenith ionospheric Faraday depth from 
day i and the other terms are the redundant compo¬ 
nent of the visibilities. Note that V/ does not provide 
contributions at high k to the power spectrum due to 
the assumption of smooth-spectrum foregrounds, and we 
are concerned only with the leakage due to Vq and Vu- 
These visibilities’ contribution to the power spectrum for 
given LST will average approximately as 


Poc |Up 


1 





(14) 









































New Limits on Polarized Power Spectra at 126 and 164 MHz 


7 


The sum may be rewritten in terms of the i = j compo¬ 
nents and the j > i component: 

=7V-h2^cos{2 ($,-$j)A2}. (15) 
i,j i>j 

In the limit where all values of 4)^ are equal, this second 
term becomes N{N — l)/2, the number of i,j pairs with 
i > j, showing that with no daily fluctuations in iono¬ 
spheric Faraday depth, there is no change in the signal. 

To estimate the level of ionospheric attenuation, we 
calculated e for the 3 LST transits described in the sec¬ 
tion above. For the observed distribution of RM at each 
of the three representative LSTs, the average attenuation 
was calculated to be e = 0.015. To obtain some measure 
of the uncertainty in this value, we generated 10,000 re¬ 
alizations of the attenuation factor for an 82 day integra¬ 
tion, drawing randomly from the RM distributions shown 
in Figure and found the average simulated value to be 
0.012 ± 0.004. Note that while we have neglected it, the 
actual spatial variation in TEC, would tend to rapidly 
increase the attenuation. The net effect is that polarized 
leakage into Stokes I is strongly decreased by this aver¬ 
aging process, as well as the observed signal in Stokes Q 
and U. 


4. RESULTS 

4.1. Features of the Power Spectra 

Figure shows the power spectra resulting from the 
above anm^ysis in I, P = Q + ill and U as a function 
of LST for Band II. Figure |7] shows the LST average 
of the power spectra of the four Stokes parameters in 
Bands I and 11. The uncertainty in these powe r spectra 
is the bootstrap error described in Section |3.6[ The ex¬ 
pected thermal noise sensitivity using the Tgyg computed 
from Section |3.5| a nd t he sensitivity calcul ations of|P ar-| 
sons et al. (|2012a|) and Pobej^et al. (20141 are shown in 
dashed, cyan lines in Figure 

We comment on some notable features of the power 
spectra. As is evident from Figure]^ the Band II Stokes 
I and V are largely devoid of features in LST, and the 
Stokes / power spectrum in Figure reproduces that 
of P14 and J14, as expected. Stokes V is very nearly 
consistent with zero, except for the slight excess below 
k ^ 0.15, which is clearly due to the one observed feature 
in Figure 1^ of unknown origin. It is also evident that 
Band I is less consistent with zero than Band II (for all 
Stokes parameters), with the excess particularly notable 
for k < 0.2. These features were noted by J14 for Stokes 
/ in the lower frequency band, and attributed to poorer 
cleaning of foreground power at lower frequencies and 
near the edge of the “wedge”. We note that a number 
of calibration leakage terms which would spill smooth- 
spectrum Stokes I into Q or U would be confined to inside 
the “wedge”, corresponding to the region blanked out in 
Figure For a f urther discussion of these effects, see 
Kohn et al. (2016). 

'there are two features in Pp worth noting. The first 
is the excess of emission at |/c||| < 0.2 hMpc”^, between 
right ascension IhOOm and 4h30m. That the excess ex¬ 
ceeds the power in the corresponding fc-bins of Pj in¬ 
dicates that it is not due to leakage from Stokes I. It 
roughly corresponds in RA with the diffuse, polarized 
power shown in B13, which is evident between RA 22h 


and 3h, and takes its minimum value at around RA 5h. 
The low RM reported in B13 would mean that this emis¬ 
sion would appear in our power spectra at low fc, as ob¬ 
served. Whether this emission should be identified with 
the B13 features would require a true imaging analysis 
to establish. 

The second feature of Figure!^ that we will comment 
on are the tracks of excess powerT)etween RA 6h and 8h, 
with |fc||| « 0.35hMpc“^. Such features could in princi¬ 
ple be generated by polarized point sources. However, 
mimicking the feature at fey « 0.345 hMpc”^ as a source 
passing through zenith (RA 6h52m Dec -30) would re¬ 
quire an intrinsic polarized flux of 100 Jy, whic h seems 
extraordinarily impl ausible, given the results of Hurley-] 


Walker et al. (2014). This survey covered 6100 square 
degrees, including the region above, and, excluding the 
brightest known sources (e.g., Pictor A, Fornax A, which 
are not good candidates) the brightest source found in 
Stokes I was 68 Jy. Similar considerations apply to the 
feature at negative k, and thus we regard all excess power 
in Q, U and V as most likely due to unidentified instru¬ 
ment systematics. 


4.2. Polarized Leakage into the EoR Power Spectrum 

As discussed in Sectionj^ the dominant form of leakage 
of polarized power into Pi in the PAPER power spectrum 
analysis comes from Q ^ I due to the primary beam 
ellipticity, as given in Equation]^ In general, calculating 
the fractional power leakage 


_ 


(16) 


depends on knowing both the primary beam and the sky. 
We have done simulations to determine this factor in 
M13, but here we also develop an approximation which 
allows us to estimate f knowing only the primary beam, 
under some assumptions about the source distribution 
on the sky. 

We approximate I and Q as Gaussian, random fields 
with mean zero (as an interferometer, we are insensi¬ 
tive to the mean value in any case). We also assume 
that (Q^) « p^{P), that is, there is an average polariza¬ 
tion fraction p of Stokes I relative to Q, where p 1 


(which is ind eed true at higher frequencies; Tucci & Tof- 
folatti]2dl2 ). Under these assumptions, we can write the 


square ot the visibility Vj, proportional to the measured 
I spectrum, as 


\Vi\^ =A+^Pi+A-^Pq, (17) 

where Pj and Pq are the true power spectra of I and 
Q, and ® denotes a convolution. The weighting factors 
A± are the contributions in power from the summed and 
differenced primary beams, defined as 


A± = J \A,, ± Ayy\^ dn. (18) 

A similar expression to Equation can be written for 
Vq, with an interchange of I and Q. Then 

^ _ p^A+ + A- A- 
A++p‘^A- A^ 


(19) 
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Fig. 6.— Power spectra P{k) from Band II, in units of mK^ Mpc“®, shown as a function of fc|| and local sidereal time (LST) in hours, 
which is the same for all panels, fey modes within the horizon for the 30 m baselines used are masked. The left panel shows Pj, the middle 
panel, Pp = Pq + Pp, and the right panel shows Py. Features in the images are discussed further in the text. 


This ratio, when multiplied by the measured Pq can be 
used as a metric to characterize the level of polarization 
leakage present in a measurement of Pj. T a.ble shows 
the v alue of this ratio for the PAPER beam (Pober et al. 
20121 for the two bands. The ratio can also be calculated 
from simulations, and on average agrees well the simpler 
estimate of Equation 

The level of leakage predicted by Equation shows 
that in the lowest k\\ bins, the I power spectrum of P14 
and J14 cannot be dominated by Q ^ leakage. The 
levels of polarized leakage are, to order of magnitude, 
10^ mK^ in Band I, and 10^ mK^ in Band II. The lack of 
a detection of polarized power (Stokes Q and U) in Fig¬ 
ure is consistent with what is known about polarized 
emission, combined wi th the attenuation factor from the 
ionosphere in Section 3.7.2 Because we are only con¬ 
sidering A:|| > 0.15 h Mpc“^, we assume the contribu¬ 
tion of low-RM diffuse emission is negligible. For point 
sources, a polarized fraction Prmiso mhz = 0.006 (lower 
than than assumed in M13) would be consistent with the 
single reported poi nt source in B13 a nd an interpreta¬ 
tion of the results of Asad et al. (2015). When combined 
with the ionosphere attenuation, this would produce a 
power spectrum below our limits, and leakage well below 
the excess present in the lowest fey bins of the Pi power 
spectra in P14 and J14. Indeed, if these levels are typ¬ 
ical, leakage would be below th e levels reported in the 
more recent PAPER-64 results of Ali et al. (2015), albeit 


under different ionospheric conditions. 


5. CONCLUSION 

We have presented the first limits on the power spec¬ 
tra of all four Stokes parameters in two frequency bands, 
centered at 126 MHz {z = 10.3) and 164 MHz (z = 7.66). 


This data comes from from a three-month observing 
campaign of a 32-antenna deployment of PAPER. These 
power spectra are processed in the same way as the un¬ 
polarized power_S£ectrum£esults have been reporte d at 


z = 7.7 (Parson s et al. 2014) and 7.5 < z < 10.5 (Ja- 
cobs et al.||20I4). We do not find a a definitive detection 


of polarized power. The limits are sufficiently low, how¬ 
ever, that the level of polarized leakage present in previ¬ 
ous PAPER measurements must be less than 100 mK^ 
at fcy ^ 0.2 /iMpc”^, below the excess found in those 
measurements. 

We additionally find that that the variation in iono¬ 
spheric RM is sufficient to attenuate the linearly polar¬ 
ized emission in these measurements by nearly a factor 
of 100. Combined with a lo wer polarized fractio n for 
point sources than assumed in Moore et al. (2013), this 
is consistent with our non-detection of polarized power. 

PAPER in the grid array configuration presented here 
is incapable of creating the high-dynamic range images 
needed to isolate polarized emission. Future work will ex¬ 
plore in more detail the effect of the ionosphere and the 
full degree of depolarization present in long EoR observ¬ 
ing seasons, as well as the full effects of polarized leakage 
from widefield beams and polarization calibration. 
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